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Introduction 

With  the  recognition  of  the  importance  of  operation  of 
military  helicopters  at  very  low  altitudes  and  very  low  velo¬ 
cities,  the  need  for  data  on  airflow  under  these  conditions 
became  apparent.  Work  was  initiated  with  Department  of  Army 
support  to  obtain  data  at  very  low  advance  ratios  typical  of 
such  low  speed  helicopter  flight.  The  original  goals  of  the 
research  were  to  try  to  obtain  a  description  of  the 
flow  field  both  time  averaged  and  instantaneous.  Because  of 
prior  experience  with  hot  wire  anemometry,  the  effort  was 
directed  towards  using  a  hot  wire  probe  to  get  all  three  com¬ 
ponents  of  velocity. 

Experimental  Method 

The  studies  of  the  rotor  wake  were  accomplished  using  a 
2.5  ft.  (.762m)  diameter  two  bladed  teetering  rotor.  The  ro¬ 
tor  was  driven  by  an  electric  motor  mounted  on  a  six  compo¬ 
nent  balance  and  placed  in  a  4  ft.  by  8  ft.  (1.22m  x  2.44m) 
wind  tunnel.  Details  of  the  installation  may  be  found  in 
appendix  A.  The  actual  studies  of  rotor  flow  began  with  the 
recognition  that  the  task  of  obtaining  and  presenting  accurate 
instantaneous  data  on  the  time  and  space  varying  velocity 
field  would  be  most  challenging.  Two  problems  presented 
themselves,  first  the  ability  of  the  hot  wire  system  to 


La 


accurately  record  the  instantaneous  velocity  components,  and 
secondly  the  large  number  of  data  points  that  would  be  needed. 
It  was  also  anticipated  that  once  the  measurement  system  was 
developed,  that  the  time-average  flow  field  data  could  be  ob¬ 
tained  either  concurrently  or  separately  from  the  non-steady 
flow  field  data. 

Probe  Frequency  Considerations 

If  one  considers  the  probe  and  hot  wire  system  first,  it 
had  been  found  that  if  the  hot  wires  were  used  to  get  large 
amounts  of  data,  considerable  wire  breakage  could  result. 

Such  breakage  would  require  extensive  recal ibration  of  re¬ 
paired  hot  wires.  After  studying  several  possibilities  a  hot 
film  probe,  TSI  1J94-60-1S,  was  chosen.  Although  the  probe 
was  considered  to  be  rugged  enough  to  ensure  long  service, 
questions  remained  as  to  the  frequency  response  characteris¬ 
tics  of  that  probe  and  the  DISA  anemometer's  used.  Since  the 
anemometers  had  a  frequency  response  of  several  kilohertz, 
attention  was  directed  to  the  hot  film  sensors.  To  determine 
the  response  needed  one  must  consider  the  frequency  content 
that  must  be  obtained  in  the  input  signal  to  the  anemometers. 

The  frequency  content  of  interest  in  the  rotor  wake  is 
that  which  would  be  seen  by  an  ideal  probe  at  a  point  in  the 
wake.  Since  the  detailed  wake  structure  consists  of  individ¬ 
ual  vorticies,  and  vortex  sheets  moving  past  the  poinx,  as 


well  as  interacting  vorticies  and  sheets  an  ideal  probe 
would  put  out  a  signal  that  would  contain  a  significant 
pulse  at  once  per  blade  passage  frequency.  With  a  two  bladed 
rotor  this  would  mean  two  pulses  per  revolution.  Since  the 
rotor  rpm  used  was  approximately  2200  rpm,  then  the  frequency 
would  be  73  hertz.  In  order  to  analyze  the  data, 
it  was  decided  that  up  to  about  twenty  harmonics  would  be 
needed,  and  hence  it  wuld  be  desireable  to  try  and  achieve 
a  useable  probe  response  out  to  over  1  kilohertz. 

Examination  of  the  methods  used  in  hot  wire  anemometry 
led  to  the  realization  that  there  were  few  if  any  useable 
means  of  determining  the  actual  frequency  response  character¬ 
istics  of  a  hot  film  probe.  Most  techniques  in  use  involved 
an  electrical  simulation  of  the  actual  problem.  Because  of 
the  need  for  good  data  at  the  high  harmonics  of  blade  pas¬ 
sage  frequency,  a  means  was  devised  to  measure  the  frequency 
response  characteristics  of  a  large  hot  film  probe  in  a  ra¬ 
ther  direct  and  simple  fashion.  This  method  is  presented  in 
appendix  B. 

Data  Handling  Considerations 

To  handle  the  data  acquisition,  reduction,  storage  and 
presentation  of  the  non-steady  data  presented  a  unique  prob¬ 
lem.  Fifteen  hundred  points  in  the  wake  were  considered 
necessary  for  a  representation  of  the  wake.  If  data  were 


taken  at  five  degree  azumuth  positions  of  the  reference 
Mane,  then  72  points  per  revolution  would  be  needed.  To 
get  a  reasonable  sample  a  minimum  of  5  readings  would  be  re¬ 
quired  at  each  azumuth  position.  Tor  the  throe  wires  this 
meant  about  l.S  million  data  points  per  test  condition.  If 
one  were  to  strive  for  greater  confidence,  then  10  or  more 
readings  would  be  required,  and  would  at  least  double  the 
data  points.  Automated  data  handling  became  a  necessity. 

Because  no  suitable  automoated  data  handling  system  was 
available  at  the  initiation  of  the  project,  a  careful  survey 
was  made  of  available  options  at  The  Ohio  State  University. 

No  suitable  large  data  processing  system  was  found  to  exist 
which  could  handle  the  mass  of  data.  Work  was  initiated  to 
utilize  an  existing  analog  to  digital  conversion  system  which 
existed  at  the  Perkins  Observatory  seme  30  miles  north  of  the 
main  campus.  Programs  were  written  and  a  tape  system  used  to 
record  test  data.  The  tape  was  then  taken  to  the  Perkins 
Observatory  for  conversion  and  then  to  main  frame  computer 
on  campus.  All  the  system  components  were  checked  and  made 
funct ional . 

The  massive  amount  of  data  being  generated  indicated 
that  essentially  millions  of  punched  cards  would  be  generated 
if  serveral  rotor  test  cases  were  mm  and  processed. 

Study  was  then  initiated  on  the  use  of  tape  oi'  floppy 


disk  data  storage,  and  on  other  data  handling  systems.  At 
this  juncture,  about  60%  through  the  program,  the  possibility 
arose  that  the  Department  of  Mechanical  Engineering  in  the 
University  would  receive  industrial  and  University  grants  to 
establish  a  computer  aided  design  laboratory  with  a  dedicated 
minicomputer  and  associated  peripherals  including  A  to  D  con¬ 
version.  Since  the  installation  of  this  new  laboratory  seemed 
immanent,  the  non-steady  flow  efforts  were  reduced  in  magni¬ 
tude  until  such  time  as  the  new  capability  became  available. 

It  was  anticipated  that  direct  coupling  of  data  output  from 
sensors  into  the  computer  could  be  accomplished  either 
through  tape  or  real  time.  Enormously  great  savings  could 
result.  Effort  on  the  study  of  rotor  non-steady  flows  was 
reduced  in  magnitude  and  work  directed  at  a  reduced  level 
towards  obtaining  time  averaged  data  on  rotor  wakes  at  low 
advance  ratios. 

Time  Averaged  Data 

The  decision  to  delay  the  instantaneous  data  part  of  the 
program  led  to  a  concurrent  decision  to  step  up  time  averaged 
flow  field  data  acquisition.  Since  the  entire  probe  position¬ 
ing  equipment,  the  force  balancing,  rotor  model,  and  hot  film 
systems  were  operational  it  was  readily  possible  to  get  aver¬ 
age  data.  Many  runs  were  made  at  advance  ratios  ranging 
from  u  =  0.02  to  0.10.  At  the  very  low  value  u  =  0.02  the 


data  were  considered  too  erratic  to  be  useful. 


from  close 


observations  it  appeared  that  the  rotor  was  operating  in 
the  "flow  breakdown"  region  reported  by  Rae  of  Washington. 
Data  were  obtained  for  various  shaft  angles  and  blade  pitch 
angles  for  both  untwisted  and  twisted  two  bladed  rotors. 

The  detailed  description  of  the  test  procedures  and  the  data 
that  were  obtained  may  be  found  in  Appendices  A  and  0.  The 
data  at  y  =  0.10  are  similiar  in  nature  to  that  reported  by 
Hevson  and  Katsoff  in  MATA  1010.  The  most  unusual  finding 
as  discussed  in  Appendix  A,  is  the  rollup  of  the  flow  that 


passes  near  the  rotor  into  two  concentrated  vortioies  sir 


lar  to  the  action  of  a  low  aspeot  ratio  wing,  with  concur 


high  axial  velocities 


At  the  conclusion  of  the  work  reported  in  Appendix  d, 
effort  was  once  again  moved  towards  the  study  of  the  non¬ 
steady  flow  field  since  a  minicomputer  dedicated  to  work 
within  the  department  became  operational.  Effort  was  then 
directed  towards  developing  a  system  which  would  couple  the 
output  from  the  hot  film  probes  directly  to  the  minicomputer 
so  that  the  entire  system  could  operate  in  an  interactive 
mode.  The  next  section  summarizes  the  work  on  the  interac¬ 


tive  svstem  which  has  been  accomplished  up  to  the  end  of  the 
third  year  and  end  of  initial  funding  of  this  study  of  rotor 
flows . 
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An  Interactive  Method  of  Measuring  Average  Instantaneous 
Velocities  of  a  Rotor  in  Very  Low  Advance  Ratios 

The  new  minicomputer  installed  in  the  department  of 
Mechanical  Engineering  at  The  Ohio  State  University  provided 
the  means  of  improving  the  system  of  average  data  acquisi¬ 
tion  and  extending  the  work  to  average  instantaneous  veloci¬ 
ties  measurement. 

The  minicomputer  is  a  DEC  PDP  11/60  and  is  well  suited 
to  the  Fortran  based  calculations.  It  has  a  250  KBYTES 
memory  (RAM)  and  64  KBYTES  are  available  to  the  user.  The 
PDP  hosts  graphics  (Tektronix  4014-1)  and  alphanumeric 
(Infoton)  terminals.  The  terminals  are  supported  by  a  high 
speed  Printronix  Printer-plotter  and  a  Tektronix  hardcopy 
unit .  Data  and  program  storage  are  provided  by  a  dual  floppy 
disc  drive  and  two  RKO  cartrige  dies  drives. 

Two  A/D  converters  are  available  for  data  acquisition: 

1)  DEC  LPA  system  which  can  sample  up  to  40  KHz  and  can  be 
multiplexed  up  to  16  channels,  with  a  multichannel  interval 
of  22  u  sec.;  2)  GENRAD  4  channel  A/D  System  (ADS)  with  a 
sampling  frequence  of  160  KHz  and  a  multichannel  interval  of 
6  u  sec. 

For  the  present  application  the  four  channel  GENRAD  ADS 
will  be  used  initially  because  of  its  higher  sampling  fre¬ 
quency  and  lower  multichannel  time  interval.  However,  if 


more  channels  may  be  needed  (for  a  fully  automated  process, 
where  the  computer  controls  the  RPM  of  the  rotor  and  the 
X-Y-Z  location  of  the  probe)  the  DEC  L?A  system  can  be  used 


and  the  22  u  sec.  error,  between  any  sequent  channels,  can 
be  comoensated  for  by  means  of  delav  lines. 

Data  Acquisition  System 

The  data  acquisition  system  will  consist  of  the  follow¬ 
ing  components  (Figure  1):  1)  Computer;  2)  A/D  converter; 

3)  Alphanumeric  terminal;  4)  Interface;  5)  Velocity  measure¬ 
ment  system  (hot-film  anemometer).  A  500  ft.  cable  (eight 
twisted  pairs  individually  shielded)  was  extended  from  the 
computer  room  on  the  second  floor  of  the  Mechanical  Engineer¬ 
ing  building  to  the  first  floor  where  the  windtunnel  is  lo¬ 
cated.  The  connection  with  the  computer  is  complete.  The  dy¬ 
namic  effect  of  the  wires  on  the  analog  signals  and  particu¬ 
larly  on  the  triggering  signal  is  yet  to  be  investigated. 

The  development  of  the  software  was  initiated.  At  this 
early  stage  of  design,  it  seems  that  two  programs,  for  aver¬ 
age  data  and  for  average  instantaneous  data  will  be  needed, 
because  of  limited  storage  space.  However  the  outline  for 
both  programs  is  similar: 

1.  As  long  as  the  data  acquisition  is  not  fully 
automated  the  process  should  be  interactive. 


2.  The  bookkeeping  should  be  done  by  the  computer, 
and  the  operator  will  be  instructed  by  the 


computer  how  to  proceed. 

3.  The  data  collection  procedure  would  be  that 
continuity  from  previous  days  will  proceed 
with  a  minimal  information  input  needed  from 
the  operator. 

4 .  The  operator  would  be  able  to  repeat  the  last 
data  point  without  disturbing  the  normal  flow 
of  operation,  to  provide  a  check. 

5.  If  more  then  the  last  point  has  to  be  repeated 
this  would  have  to  be  done  by  interrupting 
the  normal  procedure. 

Average  Data  Collection 

Mean  value  data  averaged  over  several  rotor  revolutions 
can  provide  average  velocity,  and  that  about  l*s  KHz  contains 
the  significant  frequency  content  of  the  wake.  For  a  cor¬ 
rect  reconstruction  of  the  wave  the  sampling  frequency  must 
be  at  least  2000  Hz.  An  array,  filled  with  the  digitized 
values,  having  the  dimension  (2000,3)  will  be  filled  in  one 
second,  which  corresponds  to  about  40  revolutions.  At  each 
location  in  the  wake  the  data  will  be,  averaged,  the  effec¬ 
tive  velocities  calculated,  and  transformed  to  windtunnel 
coordinates.  The  three  components  of  velocity  (Vx,  Vy,  V z) 
and  the  location  of  the  probe  (X-Y-Z)  will  be  stored  on  the 
computer’s  hard  discs.  After  the  whole  case  will  be  fin¬ 
ished  the  file  will  be  transfered  to  the  large  ADMAHL  470 
where  the  SPEAKASY  language  would  be  used  for  the  vector 
plots . 
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Figure  1.  DATA  ACQUISITION  SYSTEM. 

Instantenous  Data  Collection 

For  the  average  instantaneous  data  collection  two  methods 
will  be  developed:  1)  Computer  on-the-line  with  the  ADS; 

2)  The  whole  case  will  be  recorded  on  a  FM  analog  recorder 
(PI-6100)  and  then  reproduced  at  the  ADS.  An  interface  should 
be  designed  and  built.  It’s  basic  specifications  and  require¬ 
ments  were  provided  to  the  electronic  lab  in  the  M.E.  depart¬ 
ment,  and  it  is  expected  to  have  it  operational  at  the  begin¬ 
ning  of  October  1979. 
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Computer  Mode  Operation 


The  inputs  to  the  interface  are: 

1.  1/REV  pulse  from  a  magnetic  proximity  pickup 
indicating  ip  =  0 . 

2.  72/REV  pulses  from  a  second  magnetic  proximity 
pickup  indicating  5  degrees  increments  in  ro¬ 
tor  rotation.  The  output  from  the  interface 
will  be  one  pulse  at  each  5  degrees  increment 
for  10  revolutions  (Figure  2). 
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Figure  2. 

The  ADS  will  be  initiated  by  the  software,  but  the  digi¬ 
tizing  process  will  start  when  the  operator  will  signal,  from 
the  interface,  that  the  probe  is  at  the  correct  location. 
Then,  the  ADS  will  be  triggered  every  5  degrees  by  the 


triggering  output  from  the  interface.  After  the  digitising 
process  has  been  completed,  the  computer  will  average  the 
data  for  every  5  degrees  increment  for  10  revolutions.  The 
three  averaged  voltages  will  be  transformed  then  in  the  usual 
manner  in  wind  tunnel  coordinates.  These  velocities  and  the 
probe  location  will  be  stored,  temporarily,  on  the  computer's 
hard  discs  until  the  case  is  completed.  Then  the  data  will 
be  stored,  permanently,  on  floppy  discs  for  further  use  and 
reference , 

Tape  Mode  Operation 

At  this  mode  of  operation,  the  interface  has  a  dual  pur¬ 
pose.  When  recording  the  data,  it  will  produce  the  trigger¬ 
ing  signals  for  each  5  degrees  increment  and  will  also  oper¬ 
ate  the  FM  tape  recorder.  The  procedure  is  as  follows:  The 
operator  will  signal,  from  the  interface,  that  the  probe  is 
at  the  correct  location,  then  the  device  will  start  the  tape 
recorder  and  after  0.5  second  delay  (to  let  the  tape  get  to 
its  steady  state  velocity)  the  analog  data  will  be  recorded 
on  3  channels.  On  the  fourth  channel  the  triggering  signals 
will  be  recorded.  After  10  cycles  have  been  counted  the 
interface  will  stop  the  recording.  In  the  reproduction  pro¬ 
cess  the  three  analog  channels  will  be  connected  directly  to 
the  ADS  and  the  fourth  channel  with  the  triggering  signals 
will  be  connected  to  the  interface  whose  role  now  will  be  the 
same  as  in  the  computer  mode. 
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Data  Presentation 

The  numerical  values  of  the  velocities  and  their  loca¬ 
tion  will  be  stored  on  floppv  discs  and  could  be  displayed 
on  any  terminal  or  printed  out  by  the  PRTNTRONIX  high  speed 
printer.  However,  the  vast  amount  of  "numbers"  (1188  points 
x  (72  x  3  numbers  per  point)  =  256608  "numbers")  makes  the 
use  of  data  visualization  a  necessity. 

For  data  visualization  the  TEKTRONIX  4014-1  graphics 
terminal  will  be  used.  This  terminal  has  a  refresh  mode  of 
operation  and  animated  pictures  of  the  wake  as  the  rotor 
passes  by,  can  be  created.  In  addition,  graphs  of  velocity 
vs.  rotor  position,  or  velocity  vs.  probe  location  can  be 
obtained  by  using  the  PLOT-10  plotting  routines  available  in 
the  PDP's  library. 


Summary 


The  effort  on  the  project  resulted  in  an  operational 
rotor  drive  and  measuring  system,  a  probe  calibration  and 
positioning  system,  frequency  response  determination,  and 
acquisition  of  average  value  data.  The  low  speed  flow  fields 
found  demonstrated  a  squirting  affect  that  could  be  seen  on 
both  sides  of  the  rotor  at  very  low  advance  rations.  Two 
non-steady  data  acquisition  systems  were  considered,  one 
developed  to  the  point  of  use  but  was  found  awkward  to  use, 
and  a  second  interactive  system  which  progressed  to  the  point 
that  initial  tests  using  it  in  a  time  average  mode  illus¬ 
trated  the  interactive  capability. 

The  next  phase  of  effort  to  be  undertaken  will  be  the 
completion  of  the  mode  control  unit  which  will  take  the  coun¬ 
ter  signals  from  the  rotor  azimuth  position  to  trigger  the 
automated  acquisition  system.  When  this  is  operational,  then 
tunnel  testing  will  resume  and  both  non-steady  and  steady  da¬ 
ta  will  be  taken.  Emphasis  will  be  placed  upon  the  sequen¬ 
tial  pictorial  presentation  of  the  velocity  data  in  a  form 
that  can  be  shown  to  groups  interested  in  seeing  the  charac¬ 
ter  of  the  rotor  flow  field. 
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Abstract 

Tha  work  revealed  that  even  at  vary 
low  advance  ratios  tha  rotor  wake  tends 
to  roll  up  vary  quickly  into  two  concen¬ 
trated  trail inq  vortices  in  a  manner 
analogous  to  a  low  aspect  ratio  winn. 

List  of  Symbols 

A  3-0  probe  coordinate  direction 

B  3-D  probe  coordinate  direction 

C  3-D  probe  coordinate  direction 


U  Component  of  r  slonq  axis  C, 
e  ft/sec 

V  Forward  fliaht  sneed  of  heli¬ 

copter,  or  wind  tunnel  free 
steam  air  velocity,  ft/sec 

V  Effective  coolino  ve loci  tv  of 

“  sensor  A,  ft/sec 

V.  Effective  coolino  velocity  of 

sensor  B,  ft/sec 

V  Effective  cool inq  velocity  of 

sensor  C,  ft/sec 


B  Anemometer  output  voltaqe.  Volts 

E  Output  of  channel  A  of  the  3-D 

*  anemometer.  Volts 

B.  Output  of  channel  B  of  the  3-D 

anemometer.  Volts 

B  Output  of  channel  C  of  the  3-D 

c  anemometer.  Volts 

k  Constant  to  express  the  deviation 

from  Cosine 

Law  of  cylindrical  anemometer 
aenaors,  dimensionless 


V  »f  Effective  coolino  velocity  of 

a  sensor,  ft/sec 

V  Component  of  t'  alono  axis  X, 
ft/sec 

V  Component  o'  t'  alono  axis  Y, 

”  ft/sec 

V  Component  of  V  alono  axis  t, 

*  ft/sec 

X  wind  tunnel  coordinate  direction 

V  Wind  tunnel  coordinate  direction 


K,  Constant  of  a  cylindrical  sensor 

1  calibration  curve 

Ka  Conatant  of  a  cvlindrical  sensor 

*  calibration  curve 

R  Rotor  radius,  ft 


*  Wind  tunnel  coordinate  direction 

a  Anqle  between  direction  of  V  and 
axis  A,  degrees 

t  Angle  between  direction  of  r  and 
axis  P ,  degrees 


U  Magnitude  of  velocity  vector, 
ft/sec 


Y  Angle  between  direction  of  V  and 
axis  C,  deorees 


0  Component  of  U  along  axis  A 
•  ft /sec 

a  Component  of  0  along  axis  B 
0  ft/sec 


o  Angla  between  direction  of  U  and 
a  cylindrical  anemometer  sensor, 

deorees 

0  Rotor  anoular  velocity,  radians/ 

sac 
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Introduction 


The  knowledge  of  the  flow  field  induced 
by  the  rotor  ie  necessary  in  any  detailed 
etudy  of  the  characteristics  of  a  hell* 
copter.  The  induced  flow  is  one  of  the 
moat  important  factors  Involved  in  the 
overall  performance,  vibrations,  structural 
reliability  and  acoustic  noise  problems 
of  rotary-wind  aircraft.  Tor  example, 
knowledge  of  the  detailed  flow  characteris¬ 
tics  results  in  a  more  accurate  prediction 
of  blade  loads  and  thus  a  more  rational 
design  of  the  blade  structure.  Such 
knowledge  could  lead  to  better  flight 
performance  of  the  hellcooter  and  reduction 
of  the  rotor  induced  vibrations.  The 
knowledge  of  the  wake  geometry  would  be 
useful  in  determination  of  the  interaction 
between  rotors  and  the  interference  of  the 
rotors  with  the  fuselage  and  other  non¬ 
lifting  bodies. 

The  need  for  accurate  experimental 
wake  data  has  increased  during  the  last 
ten  years  both  because  of  the  expanded  use 
of  helicopters  and  because  of  the  develop¬ 
ment  of  computerised  numerical  methods  for 
computation  of  rotor  wake  structures 
(Reference  1) .  These  data  are  needed 
since  the  computerized  methods  used  reguire 
a  prior  knowledge  of  the  wake  geometry. 

The  data  would  be  useful  to  determine  the 
accuracy  of  ’-he  methods  developed  bv 
comparison  of  the  measured  and  the  computed 
flow  fields. 

Rotor  hover  and  forward  flight  perfor¬ 
mance  and  wake  geometrv  have  been  studied 
by  many  authors,  and  the  number  of  theore¬ 
tical  and  experimental  works  is  relatively 
large.  However,  little  is  known  about  the 
transition  flight  speed  (0<V/flR<0. 1) 
performance.  Very  large  wabe  distortion 
and  very  severe  vibratory  loads  occur  at 
these  speeds.  The  purpose  of  this  work 
was  to  develop  the  instrumentation  and  the 
experimental  methods  and  to  measure  the 
induced  velocities  in  the  vicinity  of  a 
small  helicopter  rotor  placed  in  a  wind 
tunnel  and  operated  at  transition  flight 
speeds.  The  measurements  were  performed 
by  a  three-dimensional  hot-film  anemometer. 
A  three-dimensional  traverse  system  was 
used  to  traverse  the  anemometer  probe 
throughout  the  rotor  wake.  The  experimen¬ 
tal  study  conducted  concentrated  on  advance 
ratios  of  u  «  0.1  and  below  and  the  data 
presented  are  time  averaged  values. 

Prior  forward  night  Wake  Studies 

The  rotor  wake  in  forward  flight  is 
more  complex  than  in  hover,  since  here  the 
RPepwtea  -ariatioe  of  the  bound  vortex 
strength  is  a  function  of  the  aximuthal 


position.  Tn  steady  state  forward  flight 
the  cycle  of  variation  repeats  itself  for 
each  and  every  rotor  revolution.  The 
wake  in  forward  flight  conaists  of  a  vortex 
sheet  and  a  tin  vortex,  like  ir  hover,  but 
both  components  are  more  complex.  "Ke 
vortex  sheet  may  be  represented  pictorial ly 
by  a  mesh  of  shed  and  trailing  vortices, 
as  shown  in  Figure  1  vhics  is  reproduced 
from  Reference  2.  "be  strength  of  the 
shed  vortex  elements  is  not  constant  along 
the  sran.  The  variation  corresponds  to 
the  aranwise  variation  of  the  hound  vortex 
which  existed  on  the  blade  when  the  vortex 
was  shed  due  to  chances  in  circulation  as 
the  lift  varies  around  the  azimuth.  In 
addition  to  the  shed  vortex  elements  there 
are  trailing  vortex  elements,  which  are 
deposited  to  the  wake  because  of  the  chance 
of  the  bound  vortex  strenath  distribution 
as  the  blade  chances  its  aslmuthal  position. 

The  rotor  wabe  in  forward  flight  does 
not  loch  like  a  simnle,  helical  sheet  of 
vorticitv  as  predicted  bv  the  momentum 
theory.  The  tip  vortex  is  laroelv  dis¬ 
torted  and  can  becone  unstable  and  even 
disappear.  Figure  2,  reproduced  from 
Reference  1,  shows  a  computer  simulated 
tip  vortex  traiectorv  in  forward  flight. 
Figure  3,  reproduced  from  Reference  3, 
shows  the  rolling  up  of  the  waVr  in  forward 
flight . 

Most  of  the  data  on  the  rotor  wake  in 
forward  flight  come  'rom  flow  visualisation 
studv.  Rome  authors  used  Pitot-tubes  to 
measure  the  velocitv  vectors,  (3,  4,  5,  * ) . 
Data  which  describe  the  rotor  wake  velocitv 
distribution  in  forward  flioht  are  limited 
(7,  8,  «) .  Fsoeciallv  rare  are  the  data 
on  transition  speed  forward,  flight,  o'  heli¬ 
copters,  which  are  needed  by  the  helicopter 
designers  because  o'  the  large  wake  distor¬ 
tion  and  severe  rotor  vibrations  which 
occur  at  these  speeds.  This  work  was  an 
attempt  to  develop  the  instrumentation  the 
experimental  methods  needed  to  produce 
these  data,  and  produce  the  data  at  low 
advance  ratios. 

experimental  Apparatus  and  Procedure 
The  Test  Rig 

The  tests  were  carried  out  on  a  2.5 
ft.  model,  teetering,  two-blade  rotor  in 
the  test  section  of  the  8  ft.  x  4  ft.  wind 
tunnel  of  The  department  of  Mechanical 
engineering  at  "he  Ohio  State  t’niversi  tv. 
The  details  of  the  rotor  are  shown  below. 

Details  of  the  Rotor 


dumber  of  Rlades  2 


Rotor  Radius,  ft  2 

Blade  Chord,  ft  0.177 

Rotor  Solidity  0.0894 

Root  Cutout,  %R  26.2 

Blade  Taper  Ratio  1 

Blade  Aspect  Ratio  7.1 

Blade  Twist  0 

Airfoil  Section 

NACA  0012 

r-  0.0032  slug- ft2 

The  collective  pitch  and  the  rotor 
shaft  tilt  ancle  are  variable.  The  speed 
of  the  rotor  is  variable  up  to  the  maxi¬ 
mum  of  3000  rpm.  This  results  in  tir 
speeds  up  to  400  ft/sec.  The  wind  tunnel 
air  velocity  is  variable  up  'to  the  maxi¬ 
mum  of  75  ft/sec.  Fiaure  4  shows  the 
test  section  and  the  location  of  the 
rotor  in  the  wind  tunnel. 

The  Hot-Film  Anemometry 

The  Anemometer 

The  rotor  induced  velocities  were 
measured  by  means  of  a  three-dimensional 
hot- film  probe  connected  to  three  01  PA 
55D05/102C  constant  temperature  anemometer 
circuits.  To  imorove  the  resolution  of 
the  three-dimensional  anemometer  and  to 
achieve  signal  levels  suitable  for  magne¬ 
tic  recording  the  outputs  of  the  anemome¬ 
ter  were  reduced  by  the  factor  of  4  by 
means  of  voltage  dividers  and  then  ampli¬ 
fied  hy  20  bv  means  of  three  DANA  2210  DC 
differential  amplifiers.  The  differential 
feature  of  the  DANA  amplifiers  was  also 
used  to  remove  the  DC  parts  of  the  anemo¬ 
meter  outputs  which  correspond  to  sero 
velocity  incuts.  The  result  of  this 
arrangement  was  that  the  anemometer  output 
voltages  for  input  velocities  of  o  to  80 
ft/sec  were  between  0  and  10  volts. 

Directional  Sensitivity  of  Cylindrical 
Sensors 

The  basic  hot-wire  or  hot-film 
anemometer  outnut  is  a  voltage  which  is 
related  to  the  fluid  velocity  approximate¬ 
ly  as 

E2  -  Kj,  ♦  KjUl/2  (1) 

where  u  is  the  velocity  input,  E  is  the 
voltage  output  and  K. ,  X,  are  constants 
determined  by  the  anBmomBter,  the  sensor 
and  the  fluid  properties. 

If  the  fluid  velocity  is  not  perpen¬ 
dicular  to  the  sensor  Eguation  (1)  does 
not  hold,  since  the  heat  transfer  from 
the  sensor  is  reduced.  It  is  convenient 
to  define  the  "effective  cooling  velocity' 


as  the  perpendicular  velocity  which 
produces  the  same  cooling  effect  as  the 
actual,  non-perpendicular  velocity. 

Several  relationships  have  been  proposed 
to  relate  the  actual  and  the  effective 
velocities.  The  simplest  relationship 
defines  the  effective  velocity  as  the 
component  normal  to  the  sensor  of  the 
actual  velocity.  This  correlation  neglects 
the  contribution  of  the  flow  along  the 
sensor  to  the  total  cooling  effect,  and  is 
not  accurate  for  deviations  larqer  than 
40°  from  the  normal  direction.  More 
accurate  correlations  have  been  suggested 
in  References  Id,  11,  and  12. 

We  used  the  relationship  proposed  in 
Reference  11,  because  it  matched  the 
experimental  data  accurately.  Tt  also 
enabled  us  to  obtain  a  closed  form  solu¬ 
tion  for  the  magnitude  of  the  velocity 
vector  and  its  directional  cosines. 

Figure  5  and  Eciuation  (2)  define  this 
relationship . 

V2ff  -  l’2  (sin2i  +  k2eos2«)  (2) 

k  is  a  dimensionless  constant.  Note  that 
if  the  flow  is  normal  to  the  sensor  S  » 

90°  and  Vpff  -  U. 

Data  Reduction  Method 

The  three-dimensional  probe  consists 
of  three  orthoqonal,  cylindrical  sensors, 
as  shown  in  Fiqure  6.  The  effective 
cooling  velocities  of  the  three  sensors 
are 

V2  ■  l,2(sin2o  *  k2cos2n)  (3a) 

V2  -  t’2(sin2P  ♦  k2cos2C)  (3b) 

P 

V2  -  D2(sin2y  ♦  k2cos2y)  (3c) 

c 

Adding  the  three  above  enuations  leads  to 

V2  +•  V2  ♦  V2  -  t,2{  (sin2a  ♦ 

sin2?  ♦  sin2y)  ♦  k2(cos2ct  ♦ 

cos2S  ♦  cos2y) )  (4) 

The  following  trigonometric  relationships 
exist  between  the  directional  angles  a, 

8,  Y . 

sin2a  *  sin28  ♦  sin2Y  ■  2  (5) 

2  2  2 
cos  a  ♦  cos  8  ♦  cos  y  ■  1  (6) 

Substituting  Fouatlon  (5)  and  Fouation  (6) 
into  Fquation  (4)  gives 

v?  ♦  V2  ♦  V2  -  D2 (2  ♦  k2)  (7) 

ADC 


and  solving  for  the  mean  velocity 


"\ 


/v2  +  V?  ♦  V2 
/  »  b  c 

V  2  ♦  v 


transformation  of  coordinates  was  required. 
The  two  sets  of  axes  are  shown  in  Figure  8. 
The  orientation  of  the  probe  is  such  that 
axis  C  lies  in  plane  X-7  and  olane  BOAD  is 
(8)  perpendicular  to  place  X-Z.  The  geometry 
of  the  probe  is  such  that 


After  substituting  Equation  (8)  into  <  FOC  »  35.26° 

Equation  (3)  and  some  algebraic  manipula¬ 
tions  one  can  write  <  EOD  -  54.74° 


sine  “ 


sinS 


siny 


/  (r^)(vT7^fpv) 

(9a) 


<  AOD  -  <  BOD  «  45.0° 

Using  basic  trigonometry,  we  now  can 
-yi  express  the  velocity  vector  u,  given  by 

-  '  means  of  its  components  U  .  IV.  and  V  , 

AD  C 

in  the  tunnel  system  of  coordinates  by 


(9c) 


where 


Equation  (8)  and  Eauation  (9)  define  the 
total  velocity  vector  in  the  probe  system 
of  coordinates,  provided  one  knows  by  some 
other  means  whether  the  angles  a,  B,  y  are 
smaller  or  larger  than  90°.  Since  the 
cylindrical  sensors  do  not  distinguish 
between  yaw  angle  6  and  yaw  angle  <180°  - 
4) ,  the  set  of  the  three  output  voltages 
of  the  three-dimensional  anemometer  can 
correspond  to  2J  ■  8  different  spatial 
orientations . 

Transformation  of  Velocities 

The  probe  used  was  a  TSI  1294-60-18 
hot-film  probe.  The  three  mutually  per¬ 
pendicular  sensors  are  inclined  54.74° 
to  the  probe  axis.  The  probe  was  placed 
in  the  tunnel  as  shown  in  Piaure  7. 
Experiments  have  shown  that  this  was  the 
optimal  spatial  orientation,  since  the 
average  velocity  vector  in  the  tunnel 
stayed  in  one  octant  for  most  of  the 
experiments.  The  above  mentioned  octant 
is  the  one  determined  by  outward  pointing 
vectors  along  sensors  A  and  B,  and  an 
inward  pointing  vector  along  sensor  C. 

Since  the  velocity  data  was  needed 
in  the  tunnel  system  of  coordinates,  a 


U  *  U  coso  U  has  direction  of  B  (13a) 

*  A 

Ujj  “  U  cosB  Ufa  has  direction  of  B  (13fc) 

U  ■  U  cosy  U  has  direction  of  C  (13c) 

c  c 

Calibration  of  the  Three-Dimensional  Probe 

The  calibration  of  the  three-dimen¬ 
sional  probe  consists  of  two  stages.  The 
first  is  to  obtain  the  voltage  versus 
velocity  curve  for  each  sensor,  ’’’he  second 
is  to  determine  the  yaw  ancle  sensitivity 
constant  k. 

Velocity  Calibration 

A  calibrator  similar  to  TP I  Model 
1125  was  modified  to  make  the  calibration 
of  the  3-D  probe  possible.  The  calibration 
procedure  is  similar  to  that  described  in 
references  (13,  14). 

Using  the  calibration  rio,  curves  were 
obtained  of  anemometer  output  volts  at 
various  velocities  for  all  three  orobes. 
since  the  subsequent  calculation  of  the 
velocities  was  done  by  a  digital  comouter, 
it  was  desirable  to  fit  a  polynomial  to  do 


the  calibration  curva  of  each  sensor.  The 
first  approximation  to  tha  calibration 
curve  of  a  cylindrical  sansor  is  Equation 
(1) .  To  achiava  better  accuracy  v*  used 
fourth  order  polynomials. 

Tha  polynostials  used  to  reduce  tha  data 
collected  were 

V  •  0.01546  E*  -  0.2395  E3  ♦  1.906  E2 

a  p  a  a 

-  3.03  Em  ♦  0.9  (14a) 


The  ’alue  of  V  •  0.10  was  used  in 
the  calculations  throuohout  this  work. 

Tt  is  interesting  to  note  that  the  value 
of  t.'  obtained  via  Fnuation  (9)  and  the 
values  of  n,  0  and  v  obtained  via  Fnuation 
(9)  are  only  sliahtlv  effected  by  large 
chanaes  of  k*.  The  velocity  chanoea  only 
0.24%  when  k*  changes  from  n.no  to  0.10. 

The  ancle  «  for  example,  chanaes  onlv  by 
1.56°  when  k^  chanaes  from  n.oo  to  O.in  for 
the  nominal  value  of  n  ■  inp,  which  is  a 
large  yaw  angle. 


0.01424  E*  -  0.2111  r?  ♦  1.677 
b  b 

1.65  tb  -  0.1 


"b 

(14b) 


V  - 

c 

T 


0.00100  T*  *  0.0152 
c 

1.01  Ec  -  0.5 


F3  ♦  0.410  E* 
c 

(14c) 


vhers  the  constants  were  obtained  by  a 
least  squares  fit  of  the  data. 


The  probe  was  recalibrated  periodically 
and  only  negligible  changes  were  observed. 


Directional  Sensitivity  Calibration 


Prove  Traverse  Fvstem 

The  purpose  of  the  traverse  system  is 
to  locate  the  hot-film  probe  accurately 
within  the  wind  tunnel  test  section,  while 
the  tunnel  and  the  rotor  are  running.  The 
system  is  based  on  a  X-Y  positioner  which 
was  built  for  another  proiect,  several 
years  ago.  The  X-Y  oositioner  is  capable 
of  moving  tKe  probe  support  16  inches  hori¬ 
zontally  and  10  inches  vertically  hv  means 
of  power  trains  moved  by  two  DC  motors. 

The  motion  is  controlled  from  a  control 
panel,  and  location  is  measured  bv  counting 
the  revolutions  of  the  cower  screws  using 
electromechanical  counters.  The  number  of 
counts  per  inch  of  motion  is  13,  so  that 
the  resolution  is  e.nfl  inches. 


The  purpose  of  this  calibration  was 
to  determine  the  constant  k  from  Fouation 
(2) .  The  enule  of  the  flow  to  the  wire 
was  set  to  qive  angle,  6,  between  90°  and 
0°.  The  calibrator  velocity  is  set  to  a 
constant  value  U,  and  the  output  of  the 
anemometer,  E,  is  recorded  for  several 
values  of  6.  Then,  using  the  calibration 
curve  of  the  anemometer,  the  effective 
coolinq  velocities  are  computed.  One  can 
now  rewrite  Equation  (2)  as 

■  ain2«  ♦  k2cos26 

"X  2 

and  obtain  the  value  of  k  usina  the 
Method  of  Least  Squares. 

This  calibration  waa  performed  three 
times  during  a  period  of  six  months,  for 
the  three  sensors  of  the  3-D  probe. 
Velocities  of  30,  50  and  00  ft/sec  were 
used.  The  values  of  k?  were  always  in  the 
range  of  0.090  to  0.102.  The  data  for 
S  <  10*  were  not  used  in  determination  of 
k*.  aince  the  disturbance  due  to  the 
probe  prongs  at  these  yaw  angles  was 
large  and  the  data  did  not  fit  Equation 
(2).  The  fit  for  angles  larger  than  10* 
was  considered  to  be  good,  and  the  fit 
f or  angles  smaller  than  10*  was 
considered  to  be  fair. 


To  move  the  traverse  assembly  along 
the  tunnel  in  the  flow  direction  the  device 
was  eguiDred  witv  wheels.  The  moving 
traverse  unit  is  equipped  with  *our  cam 
followers  which  slave  it  to  a  12  ft  lone 
track  made  o*  standard  steel  angle  fixed  to 
the  tunnel  floor.  The  X-V  positioner  can 
be  shifted  +  6  inches  perpendicularly  tc 
the  track  bv  relocating  the  cam  followers, 
which  is  accomplished  in  seconds,  "'•be 
motion  along  the  tunnel  is  controlled 
external  to  the  tunnel .  •'•he  vertical  range 

of  the  X-Y  positioner  was  extended  as 
necessary  by  positioning  the  probe  at  four 
different  vertical  locations  on  the  trav¬ 
erse.  h  continuous  vertical  travel  of  4^ 
inches  was  achieved.  The  soace  covered  by 
the  traverse  system  is  28  in.  x  100  in.  x 
40  in.,  and  it  can  be  increased  even  fur¬ 
ther  by  relocating  the  track.  This  range 
allows  flow  beneath  the  rotor  to  be 
studied.  Pv  resetting  the  rotor  Pitch  the 
same  traverse  system  can  be  used  to  studv 
flow  above  the  rotor. 

Velocities  were  measured  at  arid 
points  located  at  3  inch  spacinns  ver¬ 
tically,  4  inch  snacinos  across  the  rotor 
in  the  X-Y  plane,  and  at.  4  inch  soacinn 
along  the  tunnel  axis.  For  the  tests  o* 
rotor  flow  near  to  and  below  the  rotor  a 
total  of  1188  grid  points  were  available. 


Fvaluation  of  Errors 


the  nanometer  which  measure  the  Pitot 
pressure  was  estimated  to  he  from  2.5% 
at  20  ft/sec  down  to  0.25t  at  70  ft/sec 


Velocity  and  Angle  Errors 

Since  the  results  of  this  work  are 
only  as  oood  as  the  data  taken,  and  the 
methods  of  reducing  the  data,  a  brief 
discussion  of  the  maqnitude  of  the  errors 
involved  on  the  work  will  he  nresented, 

A  more  detailed  analysis  can  he  found  in 
reference  10. 

The  velocity  measurements  are  con¬ 
sidered  to  have  an  average  error  of  0.93 
ft/sec  over  the  entire  velocity  calibra¬ 
tion  range  of  10  to  70  ft/sec.  Detailed 
study  of  the  effect  of  tunnel  and  nrohe 
vibration  revealed  that  vibration  had  a 
negligible  effect  on  the  accuraee  of 
velocity  measurement.  The  estimated  error 
in  the  measurement  of  angles  with  the 
probe  indicated  that  the  maximum  errors 
were  within  +  2°  in  the  ranges  of  anales 
useful  in  this  study  'or  the  entire  range 
of  velocities  considered. 

The  errors  in  determining  the  velo¬ 
cities  in  the  x,  y,  r  coordinate  systems 
deoend  uoon  the  velocity  and  the  anole 
errors  of  the  nrohe  data.  If  the  worst 
errors  occur  simultaneously,  an  error  of 
1.9  ft/sec  could  occur.  The  averaae  error 
was  considered  to  be  l.o  ft/sec  or  less. 

Other  Sources  of  Error 

The  three  sensors  of  the  3-D  nrohe 
arc  not  at  the  same  ooint  in  the  wake, 
because  of  their  finite  size.  Tbev  'it 
within  a  0.3S  inches  diameter  sohere.  The 
error  due  to  this  arrangement  will  be 
small,  axcept  for  the  case  when  a  vortex 
core  or  a  vortex  sheet  cross  this  sphere. 

The  error  of  the  nrobe  location  in  the 
tunnel  <X,Y,Z)  is  a  result  o'  the  deflec¬ 
tion  of  the  orobe  holder  and  the  nrobe, 
imperfection  of  the  tunnel  floor  on  which 
the  traverse  mechanism  moves,  errors  of 
coordinate  measurements,  etc.  It  was 
estimated  to  be  less  than  +  0.25  inches  in 
each  coordinate  direction.- 

The  angular  error  due  to  misalignment 
of  the  probe  with  respect  to  the  tunnel 
walls  waa  estimated  to  be  ♦  0.5°.  However, 
this  error  is  included  in  the  analysis  of 
Section  3.5.2,  since  the  same  problem 
arises  when  the  nrobe  is  mounted  on  the 
velocity  calibrator.  The  angular  align¬ 
ment  waa  done  with  an  accurate  mercury 
balanced  indicator. 

The  error  of  the  tunnel  velocity  meas¬ 
urement  due  to  the  finite  resolution  of 


Changes  of  the  room  temperature  miqht 
cause  errors  in  the  tunnel  velocity  measure¬ 
ments  (Pitot  tube) ,  and  in  the  wake  velo¬ 
city  measurements  (anemometer).  However, 
we  attempted  to  minimize  this  source  of 
error  by  keening  the  temperature  constant 
within  78°  +  2°.  It  was  done  bv  controlling 
the  number  of  onen  windows  in  the  room 
where  the  tunnel  was  located. 

The  PPM  o'  the  rotor  was  steady.  It 
was  measured  hv  an  electronic  counter,  hv 
countina  the  number  o'  revolutions  over 
the  period  of  10  seconds,  mhe  accuracy  of 
this  measurement  was  better  than  1%. 

The  error  o'  the  commuted  advance 
ratio,  which  results  from  the  errors  of 
the  tunnel  velocitv  and.  the  rotor  Rpw  was 
computed  to  be  from  3.5*  at  tunnel  velocity 
of  20  ft/sec  down  to  1.3*  at  tunnel  velo¬ 
city  of  7Q  ft/sec. 

The  errors  of  the  rotor  shaft  tilt 
angle  and  of  the  collective  pitch  anqle 
were  estimated  to  be  within  the  limits  of 
+  0.5  . 


Test  Results 

The  parameters  of  the  fliaht  condi¬ 
tions  tested  in  this  seauence  of  tests 
were : 


Blade  Tip  Sneed,  ft/sec 
Advance  Ratio 

Collective  Pitch  Anole,  degrees 
Rotor  Shaft  Tile  Ancle,  degrees 
Conino  Angle,  degrees 


300 

O.lOr 

0.06: 

8 

8 

2 


0.08: 

0.04 


The  outnut  o'  each  anemometer  was  read 
as  a  time  average  using  a  v’estinohouse  PX 
161  RMS  voltmeter  and  using  a  Fluke  8000A 
digital  voltmeter. 


The  grid  used  'or  measuring  flow  in 
the  wake  below  the  rotor  consisted  of  1188 
points,  167  o'  which  were  not  used,  because 
they  were  too  close  to  the  rotor.  Eicure 
9  shows  the  grid  points  in  the  Y-Z  plane 
for  the  case  of  u  ■  0  at  the  centerline 
X/R  ■  .27.  Eighteen  such  planes  were 
equally  spaced  between  X  »  1.60  and  X  ■ 

-  2.67,  although  not  all  are  presented. 

Several  types  of  olots  were  used  to 
interpret  the  results.  Comoonents  of  U  in 
the  X-Y,  y-z  and  Z-X  planes  were  plotted. 
This  was  dene  on  a  CRT  terminal.  Downwash 
velocity  plots  were  also  made.  The  data 


that  follow  ara  averaqe  or  mean  velocity 
data  and  arc  not  instantaneous  values. 

Data  are  presented  for  advance  ratios 
of  0.04,  0.06,  0.08,  and  0.10.  Tt  should 
he  noted  that  the  rotor  Is  not  trimmed, 
hut  simply  set  at  the  collective  pitch  and 
shaft  angle. 

Figure  10  presents  the  data  for 
u  “  0.10  for  the  transverse  Y-7.  Diane, 
which  are  planes  similar  to  the  reference 
diagram  In  figure  9.  The  exnected  rollup 
as  the  flow  moves  downstream,  (a«  X/R  gets 
increasingly  negative)  can  he  observed 
from  the  data.  Figure  11,  foV  u  »  .10 
deoicts  the  various  longitudinal  X-Z  planes 
showing  the  downwash  influence  of  the 
rotor.  Figure  12,  for  u  -  o.in  and  X-v 
nlanea  shows  the  flow  field  when  viewed 
from  "above”  the  rotor  and  the  lateral 
flow  due  to  roll-un  can  he  seen. 

For  u  -  0.08  Figure  13  presents  the 
transverse  v-Z  planes.  Figure  14  presents 
the  longitudinal  X-Z  planes  and  Figure  15 
presents  the  X-Y  nlanea. 

For  u  ■  0.06  Figure  16  presents  the 
transverse  Y-Z  nlanea.  Figure  17  presents 
the  longitudinal  X-Z  planes,  and  Figure  18 
presents  the  X-Y  planes. 

For  u  0.04  Figure  1R  presents  the 
transverse  Y-7  planes.  Flours  20  presents 
the  longitudinal  X-Z  planes  and  Finure  21 
presents  the  X-Y  nlanes. 

Due  to  limitations  on  time  data  for 
the  above  the  rotor  are  not  included  in 
this  presentation.  Data  acouired  in  the 
region  up  to  0.6  Z/R  above  the  rotor  will 
be  presented  in  the  final  report  summariz¬ 
ing  this  work. 

Discussion  of  Test  Results 
Advance  Ratio  of  0.10 

If  one  examines  the  data  presented  in 
Figure  10  the  expected  vortex  roll -up  of 
the  rotor  wake  can  readily  be  seen.  The 
picture  of  the  wake  velocities  in  the 
vertical  cross  section  are  similar  to 
what  one  would  expect  from  a  low  aspect 
ratio  wing.  The  flow  field  when  viewed 
frost  above.  Figure  12,  likewise  show  an 
expected  inflow  towards  the  rotor  center 
above  the  rotor,  changing  to  an  outflow 
below  the  rotor.  The  side  view  of  the 
rotor  flow  shown  in  Fiqure  11  likewise 
appears  somewhat  as  one  would  anticipate, 
with  up  flow  in  the  plane  beyond  the 
rotor  tip  and  down  flow  under  the  rotor. 


The  results  appear  similar  to  those  of 
Heyson  and  Katsofz  (4). 

At  an  advance  ratio  of  u  ■  0.10,  tb,. 
the  picture  that  has  been  obtained  apne  w  s 
similar  to  what  had  been  expected. 

Advance  Ratio  of  0.06 

Examination  of  the  data  for  u  ■  0.06 
as  shown  in  the  vertical  planes  in  Flours 
16  reveals  a  vortex  rollup  pattern  similar 
to  what  would  expect  behind  a  low  aspect 
ratio  wino.  Tn  this  particular  case  the 
intensity  of  the  vortex  developing  aft  of 
the  advancino  side  (on  the  right  hand  side 
of  the  figure)  seems  to  be  qreater  than 
the  vortex  developed  aft  of  the  retreating 
side . 

The  flow  patterns  when  viewed  from  the 
side  of  the  rotor.  Figure  17  show  an 
anticipated  un  flow  in  the  plane  outside 
the  rotor  tip  on  the  advancing  side,  g  ■ 

-  1.33.  Just  inside  the  rotor  tin  atR80t 
radius,  a  strong  downwash  pattern  develops, 
and  also  closer  to  the  hub  27»  radius  the 
strong  downwash  pattern  can  be  seen.  When 
we  reach  the  centerline,  at  0«  radius,  the 
downwash  pattern  is  evident,  but  it  is  not 
as  strong.  At  the  opposite  side  or  re¬ 
treating  side  at  X/R  ■  0.53,  a  highly 
concentrated  region  of  downwash  is  aaain 
evident.  The  flow  beyond  the  tip  at  the 
retreating  side,  X/R  »  1.33,  is  also  up¬ 
wards  .  The  view  from  above  the  rotor 
reveals  the  anticipated  inflows  above  the 
trailina  vortex  structure,  Z/R  ■  -0.2.  In 
the  plane  of  the  rotor,  Z/R  ■  0.0,  and 
below  the  rotor  one  can  see  the  outflow  in 
the  region  of  the  trailing  vortex. 

The  advance  ratio  data  for  y  »  0.08 
and  u  ”  0.04  are  somewhat  similar  in 
appearance  to  the  u  *  0.10  and  u  -  0.06 
just  discussed. 

Interpretation 

The  data  obtained  at  u  *>  0.10  reveal 
flow  patterns  similar  to  what  was  anticipa¬ 
ted  and  has  been  seen  in  prior  studies  at 
similar  advance  ratios  and  higher  advance 
ratios  (4).  downwash  is  similar,  and 

the  vortex  patterns  shown  when  viewed  from 
the  rear  of  the  rotor  are  expected.  The 
side  view  does  not  reveal  anvthino  unex¬ 
pected  . 

The  advance  ratio  of  u  ■  0.06  and 
other  low  advance  ratios  tested  reveal  a 
difference  which  is  believed  to  be  of 
considerable  significance.  If  one  examines 
various  cuts  of  Figure  16  as  well  as  other 


vertical  cut*  through  th«  flow  In  Figure 
17,  it  becomes  clear  that  the  flow  near 
the  rotor  tend*  to  initiate  roll  ur>  into 
two  concentrated  vorticea  in  the  rotor 
plane  and  not  juat  downstream  of  the  rotor. 
Xt  appeara  that  a  algnlf leant  oortion  of 
the  total  flow  to  the  rotor  la  "encased" 
within  these  vortices  as  they  develop. 

If  we  examine  the  magnitude  of  the  longi¬ 
tudinal  component  in  Fiaure  17  as  well  as 
the  magnitude  of  the  vertical  velocity 
eoaiponents  it  can  be  observed  that  the 
aMignitudea  are  significantly  increased,  up 
to  doubled  in  aise.  If  one  recalls  from 
invlacid  flow  theory,  if  a  vortex  stretches 
out,  a  decrease  in  diameter  results  and  an 
increase  in  tangential  velocity  follows. 

Now  if  the  vortices  shown  in  Figure  16  can 
be  Interpreted  to  be  streamlines  then  the 
fluid  once  contained  within  the  rolled  up 
region  will  tend  to  stay  there.  Since 
continuity  must  be  preserved,  as  the  vortex 
rolls  up  and  diameter  contracts  the  axial 
(continuity  component)  velocity  must 
increase.  In  a  gross  pictorial  wav,  the 
flow  coming  to  the  rotor  seems  to  wrap  up 
and  through  a  "porous"  wing  (actuator 
disk)  and  while  whirling  la  aciuirted  aft 
in  concentrated  streams  almost  like  tooth¬ 
paste  squirted  out  of  a  toothoaate  tube. 

Study  of  data  at  u  »  0.04  and  0.08 
reveals  similar  behavior.  By  the  time 
U  *  0.10  is  reached  this  pronounced  squirt¬ 
ing  is  no  longer  as  evident. 

The  potential  implications  of  these 
fundings  could  be  great  if  further  data 
continue  to  verify  their  existence.  It  is 
well  known  that  large  chanqes  in  yaw  con¬ 
trol  (tall  rotor  performance)  occur  with 
small  changes  in  crossflow  or  forward 
flight  at  very  low  flight  speeds.  Like¬ 
wise  significant  changes  to  horiiontal 
stabilisers  or  stabllators  operating 
characteristics  occur  durinq  small  changes 
in  wind  velocity,  direction  or  magnitude 
when  operating  a  helicopter  at  low  flight 
speeds. 

Prom  the  data  obtained,  the  highly 
concentrated,  high  velocity  flows  generated 
could  impinge  upon  the  tail  rotor,  tail 
surfaces  and  the  tail  boom,  and  depending 
on  the  unique  angle  of  yaw,  pitch,  wind 
heading,  c.g.  position,  or  side  winds 
such  interaction  could  result  in  sudden 
and  significant  changes  to  the  forces 
generated  by  the  aft  located  aerodynamic 
surfaces  and  tall  rotors. 
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Velocity  vectors  in  Y-I  plan, 
for  u  »  0.04.  Coordinates 
depict  probe  position  only. 
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ABSTRACT 


A  method  to  measure  the  frequency  response  of  hot-wire 
and  hot-film  anemometers  directly  was  developed.  It  is  baaed 
on  the  comparison  of  the  input  (velocity)  frequency  spectrum 
to  the  anemometer  with  the  output  (voltage)  frequency  spectrum 
it  produces.  The  input  spectrum  was  determined  by  a  reference 
anemometer.  The  response  of  this  anemometer  was  shown  to  be 
flat  by  an  indirect,  electrical  test. 

Very  good  agreement  was  observed  between  the  results  of 
the  direct  and  the  indirect  tests  for  several  sensors.  Testing 
of  an  anemometer  using  the  new  method  takes  only  few  minutes 
and  it  determines  the  frequency  response  curve.  The  method  can 
handle  any  type  of  anemometers  or  sensors,  since  it  is  indepen¬ 
dent  of  the  electronic  circuitry  of  the  anemometers  and  the 
heat  transfer  mode  of  the  sensors. 

Since  the  frequency  response  curve  of  the  tested  anemo¬ 
meter  is  determined  in  a  few  minutes,  the  method  developed  is 
an  efficient  tool  for  frequency  optimization  of  anemometers. 

This  applies  mainly  for  hot-film  anemometers,  for  which  the 
square  wave  testing  method  leads  poor  results  as  compared  to 
hot-wire  anemometers. 
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1.  Introduction 

The  method  to  be  described  here  was  developed  to 
facilitate  the  measurement  of  the  three-dimensional 
velocity  vector  in  the  wake  of  a  small  helicopter  rotor 
placed  in  a  wind  tunnel.  To  prevent  mechanical  failure 
of  the  expensive  3-D  probe,  it  was  decided  to  use  large 
diameter  hot-film  sensors.  Since  the  frequency  response 
of  a  hot-film  anemometer  which  uses  large  sensors  is 
low,  some  of  the  high  harmonics  of  the  fluctuating  flow 
were  far  beyond  the  flat  range  of  the  instrument.  Accu¬ 
rate  measurement  of  these  harmonics  required  accurate 
measurement  of  the  frequency  response  of  the  anemometer. 

We  first  tried  to  use  one  of  the  known,  indirect, 
methods  of  anemometer  frequency  response  measurement. 

The  most  common,  method  is  the  electronic  square  wave 
test,  which  is  built  into  many  modern  anemometers.  This 
method  was  described  in  detail  by  Freymuth  (1967) ,  Frey- 
muth  and  Fingerson  (1977)  and  others.  It  works  well  for 
constant  temperature  hot-wire  anemometers,  determines  the 
cutoff  frequency  of  the  instrument  and  serves  as  a  tool 
for  frequency  optimization.  However,  for  hot-film 
anemometers  the  theory  for  square  wave  test  has  not  been 
developed,  and  one  cannot  use  it  as  an  accurate  tool 
neither  for  optimization  nor  for  estimation  of  the  cutoff 
frequency. 

The  second  indirect  method  available  is  the  sine  wave 


2 


test,  described  by  Weidman  and  Browand  (1975) ,  Freymuth 
(1977)  and  Freymuth  and  Fingerson  (1977) .  We  used  this 
method  to  determine  the  accuracy  of  the  direct  method 
developed.  This  method,  although  quite  simple  if  the 
techniques  suggested  in  Section  3  of  this  paper  are  used, 
is  not  practical  for  frequency  optimization. 

Frequency  optimization  is  a  trial  and  error  procedure, 
and  several  trials  are  required  until  the  optimum  is 
reached.  Since  the  sine  wave  test  is  quite  time  consum¬ 
ing,  performing  it  several  times  may  take  many  hours. 
Another  disadvantage  of  this  method  is  its  dependence  on 
the  knowledge  of  the  form  of  the  anemometer  transfer 
function.  If  the  assumed  function  is  not  accurate,  the 
results  will  be  poor. 

Because  of  the  above  disadvantages  of  the  indirect 
methods,  we  decided  to  use  a  direct  method  which  would 
be  independent  of  the  electronics  of  the  particular 
anemometer,  will  use  a  velocity  signal  as  its  input  and 
will  be  fast  enough  to  make  "on-line"  frequency  optimi¬ 
zation  of  anemometers  possible. 

Prior  to  developing  our  method  we  tried  to  adopt  one 
of  the  known  direct  methods.  However,  we  did  not  find 
one  which  would  fulfill  the  above  requirements.  Measuring 
the  response  by  placing  the  sensor  in  a  vortex  street 
shed  from  a  cylinder  (grant  1968,  Perry  and  Morrison  1971) 
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was  considered  too  complicated  and  time  consuming,  since 
a  new  combination  of  velocity  and  cylinder  diameter  is 
needed  for  every  point  of  the  frequency  response  curve. 
The  method  of  vibrating  the  sensor  on  a  high  frequency 
vibrating  table  (Perry  and  Morrison  1971)  is  limited  to 
low  frequency  and  low  perturbation  velocity  by  the  in¬ 
ertia  forces.  The  problem  of  sensor  resonance  may  also 
arise.  The  method  of  audio-frequency  modulated  micro- 
waves  directed  at  the  sensor  (Kidron  1966)  was  not 
actually  a  direct  method  since  the  input  was  not  velocity. 
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Direct  Determination  of  Hot-Wire  or  Hot-Film  Anemometer 
Frequency  Response 

The  method  developed  and  presented  here  is  based  on 
the  following  equation  known  from  signal  theory  (Doebelin 
1975).  Let  G(ioj)  be  the  sinusoidal  transfer  function  of 
a  dynamic  system,  (iw)  the  input  spectrum  to  that  system 
and  Q0  (iti))  its  output  spectrum.  Then  one  can  write: 

««•>  - 


Note  that  for  a  particular  frequency  the  three  above 
functions  are  just  complex  numbers.  If  one  is  not 
interested  in  the  phase  angles,  he  can  write: 

j  G  (io>)  |  =  l&iiaU  (2) 

!q£ (iw) | 

which  relates  the  absolute  values  (amplitudes)  of  the 
functions.  Note  that  the  plot  of  | G (ioi)  |  versus  the 
frequency  to  is  the  frequency  response  curve  of  the  dy¬ 
namic  system. 

The  logic  behind  the  method  is  simply  to  excite  the 
anemometer  with  a  signal  of  known  frequency  spectrum,  to 
measure  its  output  spectrum  and  to  obtain  the  frequency 
response  curve  by  means  of  equation  (2) . 

The  fluctuating  velocity  signal  was  obtained  using  a 
small,  radial,  high  speed  fan  taken  out  of  a  vacuum 
cleaner.  It  produced  average  velocities  up  to  25  m/s 
and  fluctuations  up  to  5  m/s  RMS.  Different  combinations 
of  the  average  velocity  and  the  RMS  value  were  available 
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at  different  locations  in  the  flow.  The  location  we 
used  is  sketched  in  Figure  1.  The  average  velocity  at 
that  point  was  15  m/s  and  the  fluctuations  were  2  m/s 
RMS.  The  frequency  spectrum  produced  by  the  fan  is 
shown  in  Figure  2.  It  was  measured  by  a  DISA  55  DO5/102C 
constant  temperature  anemometer  with  a  TS1-1210-T1. 5 
hot-wire  probe.  The  response  of  this  system  was  found  to 
be  flat  up  to  3.5  kHz  using  an  indirect ,  electric  method 
to  be  described  in  section  3.1.  The  portion  of  the  spec¬ 
trum  curve  for  frequencies  above  3.5  kHz  was  corrected 
according  to  the  frequency  response  curve  of  the 
anemometer,  and  smoothed  to  obtain  Figure  2.  Same  fre¬ 
quency  spectrum  curve  was  obtained  later  using  a  TS1  1050 
constant  temperature  anemometer  with  the  same  hot-wire 
probe.  Mote  that  the  spectrum  had  significant  content 
up  to  10  kHz. 

The  fan  was  found  to  be  a  stable  source  of  signal. 

The  average  velocity,  the  RMS  value  and  the  frequency 
spectrum  curve  did  not  change  with  time. 

The  measurement  of  the  spectra .  was  performed  by  a 
Nicolet  UA-500  Ubiquitous  Real  Time  Spectrum  Analyzer 
Calculator.  This  very  fast  and  accurate  analyzer  is 
linear  to  within  +  1/4  db  and  is  capable  of  spectrum 
averaging.  The  spectrum  curve  is  displayed  on  a  display 
scope  and  can  be  plotted  using  a  X-Y  plotter.  Accurate 
measurement  of  amplitudes  is  achieved  by  a  digital  meter 
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with  0.1  db  resolution.  Two  curves  can  be  displayed  on 
the  display  scope  at  once,  since  the  instrument  has  a 
digital  memory.  Both  linear  and  logarithmic  coordinates 
are  available. 

The  first  step  of  the  frequency  response  measurement 
procedure  was  to  obtain  the  input  spectrum  and  to  store 
the  curve  in  the  memory  of  the  spectrum  analyzer,  as 
described  above. 

The  second  step  was  to  measure  the  output  spectrum 
of  the  anemometer  being  tested.  It  was  done  by  locating 
the  probe  at  the  same  location  in  the  flow  where  the 
reference  probe  was  before.  The  accurate  location  was 
found  using  the  output  of  the  tested  system  which  should 
read  the  average  of  15  m/s.  The  output  spectrum  was  then 
analyzed  and  displayed  on  the  screen  together  with  the 
input  spectrum.  The  calculation  of  the  frequency  response 
was  simplified  using  decibel  coordinates  for  the  amplitudes. 
Equation  (2)  expressed  in  decibels  becomes: 

db  -  20  logio!  G  (iw)  |  *  20  log  ioi00  (iu>)  j -20  logio  |  (iu>)  (  (3) 

and  the  frequency  response  of  the  anemometer,  in  db,  is 
simply  the  difference  between  the  output  and  the  input 
spectra,  in  db.  Before  the  subtraction  was  performed, 
the  two  spectrum  curves  were  shifted  vertically  so  that 
their  low-frequency  regions  coincided.  (Note  that  every 
anemometer  has  flat  response  up  to  some  frequency.  If  the 
curves  appear  vertically  shifted,  it  is  because  the  static 
gains  of  the  two  sensors  are  different.) 

Figure  3  shows  the  spectra  for  the  DISA  55  D05/102C 
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anemometer  with  a  hot-wire  sensor  (sensor  2,  see  Table  1 
for  details) .  The  upper  curve  is  the  input  spectrum 
obtained  with  the  T51  1210-T1.5  probe.  The  lower  curve 
is  the  output  spectrum  of  the  tested  system.  Curve  2  in 
Figure  4  is  the  plot  of  the  difference  of  the  two  curves 
It  is  the  directly  determined  frequency  response  for 
sensor  2. 

Two  hot-wire  probes  and  two  hot-film  probes  were 
tested  using  the  above  described  method.  The  same 
anemometer  and  signal  source  were  used  for  all  the  tests 
Details  on  the  sensors  may  be  found  in  Table  1.  The 
solid  curves  in  Figure  4  are  the  directly  determined 
frequency  responses.  The  circles  are  data  obtained  by 
an  electric,  indirect  test  presented  in  section  3.  The 
agreement  between  the  two  methods  is  considered  very 
good. 

Sensors  3  and  4  were  two  identical  hot-films  on  a 
three  dimensional  probe.  They  were  connected  to  two 
identically  adjusted  channels  of  the  anemometer.  However 
the  response  of  sensor  3  was  about  2  db  higher  than  that 
of  sensor  4  for  frequencies  above  400  Hz.  We  then  used 
the  developed  method  to  get  equal  responses  from  these 
two  sensors.  Figure  5  shows  the  output  spectra  of  the 
two  systems  before  the  adjustments  were  made.  Figure  6 
shows  the  situation  after  the  overheat  of  sensor  4 
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TABLE  1 

Physical  Constants  for  the  test  sensors. 


Hot-Wires,  all  on  TSl  1210  probe 

Reference  Sensor  Sensor 


Parameter  Sensor  1  2 


Material 

Tungsten, 
Platinum  coated 

10%  Platinum 
Rhodium 

10%  Platinum 
Rhodium 

Diameter  (nun) 

0.0038 

0.010 

0.010 

Length  (mm) 

1.25 

1.37 

1.77 

Overheat 

1.30 

1.13 

1.23 

Hot-Films,  all  on  TSI  1294  probe 


Platinum  Film  on 
Material  Quartz  Substrate 

Diameter  (mm)  0.152 

Length  (mm)  2 . 0 


Overheat 


1.39 
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was  increased  from  1.39  to  1.44.  Now  the  response  of 
sensor  4  was  higher  than  the  response  of  sensor  3  by 
1  to  2  db.  Next  the  overheat  ratio  was  reduced  to 
1.42,  and  Figure  7  shows  that  the  difference  between 
the  two  sensors  was  less  than  0.5  db.  The  whole 
procedure  took  less  than  10  minutes. 


v 
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Indirect  Determination  of  Hot-Wire  and  Hot-Film  Anemometer 
Frequency  Response 


An  indirect  method  for  determination  of  the  hot-wire 


and  hot-film  anemometer  frequency  response  utilizing  the 
transfer  functions  of  Weidman  and  Browand  (1975)  was 


evaluated  and  used  extensively.  This  method  was  used 
for  two  purposes.  First,  it  was  used  to  determine  the 
response  of  the  reference  anemometer  by  means  of  which 
the  input  spectrum  was  measured.  Second  purpose  was  to 
generate  data  for  the  four  test  sensors  and  to  compare 
them  with  the  results  of  the  direct  tests. 


Hot-Wires 

Weidman  and  Browand  (1975)  have  shown  that  if  the 
reactive  elements  are  neglected,  the  frequency  response 
of  a  constant  temperature  hot-wire  anemometer  was: 


where  u  is  the  perturbation  velocity  input,  e  is  the 
voltage  output,  and  Kw,  (  are  system  constants. 

If  a  sinusoidal  voltage  efc  is  applied  across  the 
wire,  the  frequency  response  becomes: 
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Note  that  equation  (5)  differs  from  (4)  by  a  constant 
and  a  first  order  numerator  term  (l  +  iuit)  which  is  precisely 
defined  by  the  time  constant  x.  Therefore,  if  one  can 
determine  the  response  described  by  equation  (5)  and 
determine  x,  equation  (4)  can  be  obtained  by  simple 
division.  If  the  frequency  responses  are  in  db  and  are 

normalized  to  their  flat,  low  frequency  portion,  the 

* 

constants  Kw  and  Kw  are  not  important. 

The  curve  in  Figure  8  is  a  db  plot  of  equation  (5) 

for  the  reference  sensor.  It  was  obtained  by  applying 

electrical  white  noise  from  a  Hewlett-Packard  HOI-3722  A 

Noise  Generator  across  the  wire  and  measuring  the  output 

spectrum.  (The  same  results  could  be  obtained  using  a 

sine  wave  generator  and  a  scope) .  During  the  test  the 

wire  was  exposed  to  a  low  turbulence  flow  of  15  m/s. 

The  dashed  line  in  Figure  8  has  the  slope  of  20  db/decade, 

characteristic  of  a  first  order  term.  Note  that  the 

accuracy  is  very  good.  The  effect  of  the  second  order 

denominator  term  on  the  slope  of  the  curve  is  negligible, 

because  w  >>  ‘/rand  the  damping  is  nearly  optimal  (C—0.7). 
n 

The  curve  obtained  by  subtracting  the  db  plot  of  a 
first  order  system  which  has  the  above  described 
asymptote  from  the  experimental  curve  of  Figure  8  is  the 
indirectly  determined  frequency  response.  It  is  shown  as 
the  dashed  curve  in  Figure  4. 
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3.2  Hot-Films 


Weidman  and  Broward  (1975)  have  also  shown  that  for 
hot-film  sensors  the  following  two  equations  apply: 


(iw)  = 


|l+  5  x  kel(x)J  ^1+iw/i  j+A-iuB 


—  (iu>)  - 

et 


Kf  *  [l+0  x  kel  (x)] _ 

1+0  x  kel(x)l  fl+iw/A  +A-iwB 


where 


kel(x)  = 


ber,*  (x)  +ibei'  (x) 
_ o  _ o _ 

bern  (x)  +ibei~  (x) 


(ks) 


Here  a  is  the  film  substrate  radius,  Ks  is  the  thermal 
diffusivity  of  the  substrate  and  Kf,  K^f  0,  H,  A,  B  are 
system  constants.  The  Kelvin  (Thomson)  functions  ber0 
and  bei0  and  their  first  derivatives  are  extensively 
tabulated  by  Nosova  (1961)  ,  and  are  also  available  in 
the  Fortran  Libraries  of  numerous  computer  centers. 

The  procedure  for  determining  the  frequency  response 
for  hot-films  is  similar  to  the  one  used  for  hot-wires. 
The  only  difference  is  that  the  numerator  term  here  is 


I l+0xkel (x)  .  A  family  of  these  curves  for  a=  0.0762  mm 
-*2 

and  Kg  —  0.84  mm  /s  is  plotted  xn  Figure  3*  For  better 
accuracy  the  indirect  test  was  performed  using  a  Weston- 

Boonshaft  and  Fuchs  711A  Frequency  Response  Analyzer 

which  measures  both  the  amplitudes  ratio  and  the  phase 


angle.  As  in  the  case  of  hot-wires,  the  accuracy  was 
very  good,  and  the  0-15  curve  from  Figure  9,  when 
subtracted  from  the  experimental  curves,  leads  to  the 
data  represented  by  circles  on  curves  3  and  4  in  Figure  4. 
The  value  of  $  was  determined  by  comparing  the  low 
frequency  portions  of  the  experimental  curves  with  the 
low  frequency  portions  of  the  curves  from  Figure  9. 

The  0*15  curves  matched  the  experimental  curves  very 
closely  up  to  80  Hz,  for  both  amplitude  ratios  and  phase 
angles. 
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4 .  Conclusions 

The  developed  method  was  found  to  be  a  useful  tool 
for  determining  the  frequency  response  characteristics 
of  anemometers.  Once  the  experimental  set  up  was  ready, 
determination  of  an  amplitude  ratio  curve  for  an 
anemometer  was  accomplished  in  few  minutes. 

Very  good  agreement  between  the  direct  tests  of 
section  2  and  the  indirect  tests  of  section  3  was 
observed,  as  shown  in  Figure  4.  The  solid  curves  (direct 
test)  and  the  circles  (indirect  test)  almost  coincide. 
Thus  the  method  developed  gives  strong  substantation  for 
the  validity  for  the  analytical  method  of  Weidman  and 
Browand. 

The  direct  method  can  be  extended  to  higher 
frequencies,  if  an  appropriate  source  of  signal  is  used. 
Our  source  had  significant  frequency  content  up  to 
10kHz .  (We  were  not  interested  in  frequencies  beyond 
2KHz ) . 

Although  we  used  a  rather  expensive  spectrum 
analyzer,  a  much  simpler  instrument  can  be  used  for  the 
direct  test.  A  general  purpose  analog  spectrum  analyzer 
used  in  the  fields  of  sound  or  vibrations  would  be  ade¬ 
quate,  however,  the  procedure  would  be  more  time  con¬ 
suming  (but  still  faster  compared  to  other  methods) . 

A  particularly  important  application  of  the  method 
presented  here  is  for  the  measurement  of  the  frequency 
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response  of  hot-film  anemometers  with  sensors  of  complex 
shape.  Even  if  the  heat  transfer  model  for  the  sensor 
is  not  known,  one  still  can  determine  the  frequency 
response  curve,  and  use  the  sensor  for  measurements  of 
fluctuating  velocities.  The  experimentally  determined 
response  curve  can  also  help  in  determining  the  heat 
transfer  model  for  the  sensor. 

We  also  found  that  the  indirect  method  described  in 

section  3  was  simple,  accurate  and  required  only  a  sine 

wave  generator  and  a  oscilloscope.  When  the  technique 

of  subtracting  db  curves  is  used ,  no  knowledge  of  the 

system  parameters  is  required  (except  a  and  K  for  the 

8 

hot-film  case) .  More  time  and  effort,  however,  are 
required  when  compared  to  the  direct  method. 
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APPENDIX 


A.  A  Different  Look  at  the  Direct  Method. 

We  have  used  the  system  dynamics  approach  throughout  the 
paper.  For  those  unfamiliar  with  this  approach,  the  following 
may  be  useful  in  understanding  the  direct  method. 

Let  the  anemometer  tested  be  represented  by  G(iu>),  the 
Fourier  transform  of  its  impulse  response  (its  transfer  function) 
and  the  reference  anemometer  be  similarly  represented  by  R(iw) . 
Let  the  input  signal  to  both  anemometers  be  Qi(iw),  the  output 
of  G  be  QG(iu>)  and  the  output  of  R  be  QR(ito)  . 

Since  G,  R,  Q^,  QG,  and  QR  are  Fourier  transforms,  one 
can  write : 


QG(iw) 

■  Q(iw) 

Qi(ioj) 

and 

QR(iui) 

-  R(iw) 

Qi (iw) 

from  which 

follows  =» 

QG(iw) 

Qu»)  ■ 

1  o^TST 

R(i  w) 

If  R(iw)  is  flat  up  to  the  highest  frequency  used,  R(iw)  ■  R^ 

*  constant  and 

x  _  „  QG(iu>)  _  „  SYSTEM'S  OUTPUT  SPECTRUM 
'  “  1  QR(iio)  "  1  REFERENCE  OUTPUT  SPECTRUM 

The  presence  of  the  frequency  independent  constant  R^  explains 
why  the  spectrum  curves  appear  shifted  on  the  screen. 

B.  Order  of  the  Anemometer  Transfer  Function 


Since  we  were  not  interested  in  frequencies  above  2KHz,  a 
second  order  denominator  term  was  adequate  to  model  the  anemo¬ 
meter.  (See  equations  4, 5, 6, 7)  For  higher  frequencies  a  third 
order  model  is  necessary,  as  explained  by  Freymuth  (1977) .  How¬ 
ever,  neither  the  direct  method  nor  the  indirect  method  would 
change  if  the  third  order  model  was  assumed,  since  they  both 
are  independent  of  the  denominator  of  the  transfer  function. 

(The  direct  method  is  actually  totally  independent  of  the  model.) 


Legends  for  Illustrations 


Relative  position  of  the  radial  fan  and  the 
hot-wire  probe. 

The  frequency  spectrum  produced  by  the  radial 
fan.  (Smoothed  Curve) 

Frequency  spectra  for  sensor  2.  The  upper  curve  ‘ 
is  the  input  spectrum  and  the  lower  the  output 
spectrum  of  sensor  2. 


The  results  of  the  direct  and  the  indirect  tests 
of  frequency  response.  The  dashed  curve  is  the 
indirectly  determined  response  of  the  reference 
sensor.  Solid  curves  1-4  are  the  results  of  the 
direct  tests  for  sensors  1-4.  The  circles  are 
the  results  of  the  indirect  tests  for  sensors  1-4 


Output  spectra  of  sensor  3  (overheat  1.39)  and 
sensor  4  (overheat  1.39). 


Output  spectra  of  sensor  3  (overheat  1.39)  and 
sensor  4  (overheat  1.44). 


Output  spectra  of  sensor  3  (overheat  1.39)  and 
sensor  4  (overheat  1.42). 


Response  of  the  anemometer  with  the  reference 
sensor  to  electrical  input  (solid  line) ;  20  db/ 
decade  asymptote  (dashed  line) . 


Plot  of  [l+6xkel(x)]  for  a*0.0762  mm  and 
1^-0. 84  mmVs. 


25  50  100  200  500  lk  7k 5 

Frequency  (Hz) 

Figure  3:  Frequency  spectra  for  sensor  2.  The  upper  curve  is  the  input  spectrum  and 
the  lower  the  output  spectrum  of  sensor  2. 
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LIST  OF  SYMBOLS 

b  Number  of  blades ,  dimensionless 
c  Blade  chord,  inch 

i  Shaft  angle,  degrees 

2 

I0  Flapping  Moment  of  Inertia,  slugs  -  feet 
2 

k  Constant  to  express  the  deviation  from  Cosine-Law  of 
cylindrical  anemometer  sensors,  dimensionless 

Veff  E^^ecti-ve  cooling  velocity  of  a  sensor,  feet/second 

V  Magnitude  of  velocity  vector,  feet/second 

X  Wind  tunnel  coordinate  direction 

Y  Wind  tunnel  coordinate  direction 

Z  Wind  tunnel  coordinate  direction 

v  Advance  ratio,  dimensionless 

0  Collective  angle,  degree 

6  Angle  between  direction  of  V  and  cylindrical  sensor 


ABSTRACT 

Tests  were  conducted  on  a  two  bladed  teetering  model  helicop¬ 
ter  rotor  and  on  a  two  bleded  rotor  with  twisted  blades  at  very  low 
advance  ratios.  Data  on  the  time-average  structure  of  the  rotor 
wake  were  taken  using  a  three  wire  hot-film  probe.  The  findings 
present  data  which  may  be  useful  in  further  understanding  the  nature 
of  the  wake  at  low  speed  flight  of  helicopters  and  may  prove  useful 
In  the  development  of  models  for  the  rotor  wake. 


1 


EFFECT  OF  CHANGING  ROTOR  PARAMETERS  ON  ROTOR 
WAKE  VELOCITIES  AT  VERY  LOW  ADVANCE  RATIOS 

H.  R.  Velkoff,  H.  Terkel,  F.  Shaio 
The  Ohio  State  University 
Columbus,  Ohio  U.S.A. 


1.  INTRODUCTION 

The  need  for  data  or  helicopter  rotor  wakes  at  low  advance 
ratios  has  increased  because  the  recent  emphasis  on  low  speed  nap- 
of-the  earth  flight.  To  provide  data  for  this  region  of  flight, 
tests  using  a  model  rotor  in  a  wind  tunnel  have  been  conducted. 

The  initial  data  obtained  utilized  a  two  bladed  model  rotor  at  ad¬ 
vance  ratios  of  0.04  through  0.10  and  were  reported  in  a  recent 
paper  (1).  In  those  Initial  tests  data  were  provided  at  only  one 
value  pitch  angle,  0  ■  8°,  and  shaft  tilt,  i  a  8°.  The  results  of 
those  initial  tests  revealed  that  the  overall  rotor  flow  tended  to 
roll-up  into  two  descreet  vortices  similar  to  a  low  aspect  ration 
wing.  A  particularly  unique  flow  characteristic  was  also  observed 
at  the  advance  ratios  of  y  ■  0.08  and  below.  This  characteristic 
was  the  pronounced  tendency  for  the  flow  to  not  only  form  the  two 
vortices,  but  to  also  present  concentrated  regions  of  inhanced 
magnitude  of  the  axial  component  of  velocity  in  the  vortex.  .Such 
localized  regions  of  high  dynamic  pressure  associated  with  the 
vortices  could  possibly  affect  the  design  of  empenage  for  helicop¬ 
ters. 


Since  data  in  the  low  advance  ratio  region  are  limited,  pri¬ 
marily  that  of  reference  1,  Heyson  and  Katsoff  (2),  and  Bowden  and 
Shockey  (3),  further  tests  were  planned  and  undertaken  to  secure  a 
wider  range  of  data.  The  data  presented  here  includes  tests  of  a 
two  bladed  rotor  at  a  given  pitch  angle  and  three  tilt  angles.  A 
twisted  blade  rotor  was  also  tested  to  evaluate  the  effect  of  blade 
twist  on  rotor  wake  velocity. 

2.  EXPERIMENTAL  ARRANGEMENT 


The  tests  were  conducted  In  a  wind  tunnel  with  a  8  foot  by 
4  foot  test  section.  The  flow  turbulence  Is  reduced  by  means  of  a 
2  to  1  contraction  ratio,  nested  1/8  inch  diameter  straws,  and  three 
layers  of  screening.  The  tunnel  with  the  turbulence  reduction  sys¬ 
tem  In  place  provides  air  velocities  up  to  75  feet  per  second. 
Further  details  on  the  wind  tunnel  can  be  found  in  reference  4. 

3.  ROTOR  CHARACTERISTICS 


The  rotors  which  were  tested  are  Indicated  below. 


Two  Bladed  Rotor 


Teetering  Rotor 

Rotor  radius 
Blade  chord 
Rotor  solidity 
Root  cutout 
Blade  taper  ratio 
Coning  angle 
Blade  aspect  ratio 
Blade  twist 
Airfoil  section  NACA 


1.25,  feet 
2.125,  inch 
bc/irR  =  0.0902 
*R  *  11.7 
■  1 
*  0° 

£  «  7.06 
c 

*  0° 

0012 

»  0.005583  slug-feet2 


Two  Bladed  Rotor  With  Twisted  Blades 


Teetering  Rotor 


Rotor  radius 
Blade  chord 
Rotor  solidity 
Root  cutout 
Blade  taper  ratio 
Coning  angle 
Blade  twist 
Airfoil  section  NACA 


1.25,  feet 
2.1875,  inch 
bc/irR  *  0.0928 
«R  «  12.08 
=  1 
-  0° 

-  8°  (From  12.08SR  to  100*R) 
001 2 

■  0.004260  slug-feet2 


The  rotors  are  driven  by  a  variable  speed  electric  motor. 
The  collective  pitch  and  the  shaft  angle  are  adjustable.  The  ro¬ 
tors  are  not  trinmed  but  are  allowed  to  flap  freely.  There  is  no 
rotor  moment  developed  in  either  of  the  rotors  tested.  Figure  1 
Illustrates  the  general  location  of  the  rotor  in  the  cross  section 
of  the  wind  tunnel. 


4.  FLOW  MEASUREMENT 


The  velocity  data  obtained  are  time  average  data  taken  with 
a  three-dimensional  hot-film  probe  connected  to  three  DISA  55005/ 
102C  constant  temperature  anemometer  circuits.  Details  of  the  hot 
wire  circuitry  used  and  the  method  of  data  reduction  may  be  found 
in  reference  1.  In  summary  each  hot  film  probe  is  calibrated  for 
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the  effects  of  changes  in  direction.  Of  specific  importance  to  ob¬ 
taining  data  at  low  tunnel  speeds  associated  with  low  advance  ratios 
Is  the  fact  that  the  hot-films  must  operate  at  very  low  probe  Rey¬ 
nolds  numbers.  At  such  conditions  the  value  of  IT  in  the  direc¬ 
tional  correction  equation 

Veff2/V2  «  sin2<S  +  k2  cos26 

2 

can  vary  significantly.  Variation  in  k  can  lead  to  errors  in  the 
computed  spatial  angles  of  the  velocity  vectors.  In  these  tests 
calibration  of  the  effect  of  low  Reynolds  number  on  k^  was  accom¬ 
plished.  A  typical  result  for  a  TSI  1294-60-18  0.006"  diameter  hot- 
film  probe  is  shown  in  Figure  2. 

5.  TRAVERSE  MECHANISM 


The  hot-film  probe  was  positioned  at  discreet  points  through¬ 
out  the  region  of  the  rotor  wake,  both  below  and  above  the  rotor. 

The  region  traversed  was  from  0.4  radius  above  the  rotor  to  0.8 
radius  below  the  rotor,  laterally  out  to  1.2  radius,  and  28  inches 
forward  of  the  rotor  centerline  and  40  inches  aft  of  the  rotor  cen¬ 
ter.  The  data  points  were  located  at  3  inches  along  rotor  axis, 

4  Inches  laterally  and  4  inches  along  the  tunnel  centerline  or  tun¬ 
nel  flow  direction.  A  total  of  1584  measurement-locations  result. 
However  the  number  of  points  is  reduced  since  data  can  not  be  taken 
in  the  vicinity  of  the  rotor  plane.  Further  details  on  the  travers¬ 
ing  system  may  be  found  In  reference  1. 

6.  ERROR  ANALYSIS 


An  extensive  discussion  of  the  errors  in  velocity  and  angle 
measurements  and  probe  location  is  also  presented  in  reference  1. 
In  general  the  error  in  angular  measurement  is  believed  to  be  of 
the  order  of  2°,  velocity  in  the  order  of  +  1  feet  per  second,  tun 
nel  speed  of  the  order  of  2%,  rotor  speed  to  within  2S,  and  probe 
position  to  within  +  0.25  inches. 

7.  TEST  RESULTS 


The  test  conditions  for  the  two  bladed  teetering  rotor  were 


Blade  tip  speed,  feet/second  300 
Advance  ratio  0.06 
Collective  pitch  angle,  degrees  8 
Rotor  shaft  tilt,  degrees  2,4,8 
Coning  angle,  degrees  0 


An  additional  test  was  run  for  a  two  bladed  teetering  rotor 
with  twisted  blades.  The  test  conditions  for  this  rotor  were 
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Blade  tip  speed,  feet/second  300 

Advance  ratio  0.06 

Collective  pitch  angle  at  75*R, 

degrees  8 

Twist  angle,  (From  12.08.R  to  100-'.")  8 

degrees 

Rotor  shaft  tilt,  degrees  8 

Coning  angle,  degrees  0 


The  data  are  presented  in  terms  of  computer  generated  vector 
plots.  Components  of  the  total  velocity,  V,  are  shown  in  X-Y,  Y-Z, 
and  Z-X  planes.  The  locations  of  these  planes  (X-Y,  Y-Z,  etc.)  is 
depicted  in  Figure  3.  It  should  be  remembered  that  the  data  pre¬ 
sented  are  average  velocity  data,  averaged  at  each  point  over  many 
revolutions  of  the  rotor,  and  that  the  rotors  are  not  trimmed,  and 
do  not  generate  rotor  moment. 

Data  are  presented  for  advance  ratio  of  0.06,  pitch  angle  of 
8  degrees,  and  shaft  angles  of  2,  4,  8  degrees.  For  the  shaft  angle 
of  8  degrees,  an  additional  condition  is  presented  with  twisted 
blades.  For  the  last  case  the  twist  is  two  degree  for  the  blade's 
length  and  the  pitch  angle  was  measured  at  751-  of  the  radius. 

Figure  4  presents  the  data  for  i  *  8°  for  the  transverse 
Y-Z  plane.  The  expected  rol-lup  as  the  flow  moves  downstream  can  be 
observed  from  the  data.  Figure  5  for  i  *  8°,  depicts  the  various 
longitudinal  X-Z  planes  showing  the  downwash  influence  of  the  rotor. 
Figure  6  for  i  *  8°  and  X-Y  planes  show  the  flow  field  when  viewed 
from  "above"  the  rotor  and  the  lateral  flow  due  to  rollup  can  be 
seen. 


For  i  «  4°  Figure  7  presents  the  transverse  Y-Z  planes. 
Figure  8  presents  the  longitudinal  X-Z  planes  and  Figure  9  presents 
the  X-Y  planes. 

For  i  «  2°  Figure  10  presents  the  transverse  Y-Z  planes. 
Figure  11  presents  the  longitudinal  X-Z  planes  and  Figure  12  pre¬ 
sents  the  X-Y  planes. 

Figure  13  presents  the  transverse  Y-Z  planes  for  i  *  S®  and 
twisted  blades.  Figure  14  presents  the  longitudinal  X-Z  planes  and 
Figure  15  presents  the  X-Y  planes. 

8.  DISCUSSION  OF  TEST  RESULTS 

Examination  of  the  data  for  u  *  0.06  was  shown  in  the  Y-Z 
planes  in  Figures  4,  7,  10,  13  reveals  a  vortex  rollup  similar  to 
what  would  be  expected  behind  a  low  aspect  ratio  wing.  In  this 
particular  case  the  intensity  of  the  vortex  developing  aft  of  the 
advancing  side  (on  the  right  hand  side  of  the  figure)  seems  to  be 
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greater  than  the  vortex  developed  aft  of  the  retreating  side.  Fig¬ 
ures  4,  7,  10  which  depict  shaft  angles  of  8,  4,  2  respectively, 
do  not  show  a  distinct  change  in  the  flow  pattern,  and  one  would 
have  to  use  a  more  thorough  statistical  analysis  to  interpret  the 
influence  of  shaft  angle  on  the  rotor's  wake.  Figure  13  reveals 
that  the  vortex  intensity  is  much  greater  for  the  twisted  blades 
but  the  general  shape  of  two  vortices  rolling  up  remains  the  same. 

The  flow  pattern  when  viewed  from  the  side  of  the  rotor 
(X-Z  planes)  Figures  5,  8,  11,  14  shows  an  anticipated  up  flow  in 
the  plane  outside  the  rotor  tip  on  the  advancing  side  (Y/R  *  1.33). 
Just  inside  the  rotor  tip  at  Y/R  *  -0.8  a  strong  downwash  pattern 
develops.  Closer  to  the  hub  at  Y/R  *  -0.267  the  downwash  pattern 
is  still  evident.  At  Y/R  *  0.0  the  downwash  pattern  is  still  ap¬ 
parent  but  reduced  in  intensity.  At  the  retreating  side  at  Y/R  * 
0.53  a  highly  concentrated  region  of  downwash  is  again  evident. 

The  flow  beyond  the  tip  at  Y/R  =■  1.33  is  upwards  similar  to  the 
plane  at  Y/R  ■  -1.33.  As  with  the  Y-Z  planes,  changes  in  shaft 
angle  do  not  reveal  a  significant  alteration  in  the  flow  pattern 
(Figures  5,  8,  11).  However,  the  twisted  blade  creates  a  larger 
longitudinal  component  (X  direction)  as  well  as  a  stronger  downwash 
(Z  direction). 

The  "squirt  effect"  discussed  in  reference  (1)  is  revealed 
in  these  experiments  again.  The  X-Z  planes  for  the  twisted  blades 
at  Y/R  *  -0:267  and- Y/R  •«  -0.533  are  actually  revealing  a  funnel 
like  behaviour.  The  flow  coming  out  of  the  rotor  between  X/R  » 

-0.8  and  X/R  ■  -1.07  is  squirted  aft  is  concentrated  streams. 

The  view  from  above  the  rotor,  Figures  6,  9,  12,  15  reveals 
the  anticipated  inflows  above  the  trailing  vortex  structure,  Z/R  s 
0.2.  At  Z/R  ■  0.0  and  below  the  rotor  one  can  see  the  outflow  in 
the  region  of  the  trailing  vortex.  For  these  planes,  again  the 
change  in  shaft  angle  does  not  show  a  significant  effect  on  the 
flow  pattern,  but  the  rotor  with  the  twisted  blades  creates  a  more 
distinct  flow  pattern  by  having  bigger  velocity  components  in  the 
x  and  y  directons. 
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FIGURE  2.  Changes  of  k  as  a  Function  of  Veff 


FIGURE  3. 


The  Spatial  Orientation  of  the 
Hot  Film  inthe  Wind  Tunnel. 


Figure  4. 

Velocity  vectors  in  Y-Z  plane 
for  V«0.06,  0-8®,  and  i-8® . 
Note  that  coordinates  depict 
probe  position  only. 
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Figure  6. 


Velocity  vectors  in  X-V  plane 
for  li«0.06,  0<*8J,  end  i-8°. 
Note  that  coordinates  depict 
probe  position  only. 
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Figure  12. 

Velocity  vectors  in  X-Y  plane 
for  y-0.06,  0"8°»  and  i-2°. 
Note  that  coordinates  depict 
probe  position  only. 
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Figure  13. 


Velocity  vector*  in  Y-2  plane 
for  y-0.06,  e?5%_  -  8#,  and 
i-8* .  Not.  tKat  coordinates 
depict  probe  position  only. 
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Velocity  vectors  in  X-Y  plane 
for  U-0.06,  8  »  -  8°,  and 

i-8°.  Note  tnat  coordinates 
depict  probe  position  only. 


